Non-universality of chaotic classical dynamics: implications for quantum chaos 
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It might be anticipated that there is statistical universality in the long-time classical dynamics of 
chaotic systems, corresponding to the universal correspondence of their quantum spectral statistics 
with random matrix models. We argue that no such universality exists. We discuss various statistical 
properties of long period orbits: the distribution of the phase-space density of periodic orbits of fixed 
length and a correlation function of periodic-orbit actions, corresponding to the universal quantum 
spectral two-point correlation function. We show that bifurcations are a mechanism for correlations 
of periodic-orbit actions. They lead to a result which is non-universal, and which in general may 
not be an analytic function of the action difference. 
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It has been appreciated for many years that sufficiently 
complex quantum systems exhibit a high degree of uni- 
versality Q : many statistical properties of their spectra 
usually fall into one of three classes, exemplified by the 
three Gaussian random matrix ensembles introduced by 
Dyson ||. It is natural to ask whether an analogous de- 
gree of universality exists in classical dynamics, and if 
it exists whether it underlies the universality observed 
in the behaviour of quantum systems. This paper sug- 
gests what the appropriate classical analog of quantum 
spectral universality should be, and gives arguments sup- 
porting the view that there is, in general, no classical 
universality underlying that of quantum systems. 

The universality exhibited by spectral properties is 
confined to statistics which are sufficiently short ranged 
in energy. This implies that the universal features are as- 
sociated with dynamics over long time scales: they may 
be associated with universal properties of the long-time 
classical dynamics, or they may be purely quantum. It is 
natural to anticipate that universal behaviour, if present 
at all, will only be manifest in properties which charac- 
terise small regions of phase space, since the large scale 
structure of phase space can clearly be non-universal. 

An example of a property which could show univer- 
sality is the statistical characterisation of the distribu- 
tion of points periodic under the N th iterate of a chaotic 
area preserving maps of the form a n+ \ = A4(a n ), with 
a = (x,p). These points are illustrated in Fig. [ljfor the 
case of the 7 th iterate of the standard map 

x n+1 = x n + p n , (1) 
Pn+l = Pn + Ksm(x n+ i) 

with K = 6 and the 11 th iterate of a modified cat map 



x n +i 

Pn+l 



1 1 

1 2 



K 
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with u n — sin(a; n + 2p n ) and K = 2. The patterns dis- 
played in Fig. [I] are complex and show very wide fluctua- 
tions in the density of points. They are clearly different, 



but local statistical properties of their fine-scale structure 
might be equivalent after scaling the coordinates to give 
the same mean density of points. 

Gutzwiller's trace formula j| gives a relation between 
classical periodic orbits and quantum spectra. The for- 
mula is exact for a small number of special systems, but 
it cannot be exact in general because it contains no infor- 
mation about the choice of quantisation procedure || . In 
[HQ Gutzwiller's formula is combined with the observa- 
tion that spectral correlations are described by random 
matrix theory to infer that the actions of long period or- 
bits are correlated. The correlations are a function of 
the action difference which was predicted to be universal 
within each of Dyson's symmetry classes, and an ana- 
lytic form was quoted for the GUE ensemble. In systems 
where Gutzwiller's trace formula is exact (for a discus- 
sion of the three known examples see [^|Jl2| ) such classcial 
correlations must certainly exist. However, in general the 
trace formula is not exact and it is thus necessary to find 
an entirely classical mechanism for such correlations. Up 
to now no such mechanism has been found. 

The paper is organized as follows. We will first ar- 
gue that the long-time, local structure of phase-space of 
chaotic systems is non-universal, using theoretical argu- 
ments and numerical experiments. We then character- 
ize the remarkably strong fluctuations of the density of 
periodic points in Fig. |lj and provide theoretical argu- 
ments as well as numerical evidence that these fluctua- 
tions are log-normal. Finally we propose an entirely clas- 
sical mechanism for periodic-orbit correlations, based on 
the statistical properties of bifurcations in the long-time 
dynamics. This mechanism gives rise to a non-universal 
and possibly non-analytic contribution to the correlation 
function. 

We will first show that the long-time, local structure 
of phase space is not universal. Our argument is based 
upon considering a particular statistic: we will consider 
the proportion P(N) of trajectories which are elliptic 
upon N iterations of area-preserving maps such as (1), 
(2). The results generalise directly to continuous time 
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systems. The stability of a point periodic under N itera- 
tions of M, a — M N (a), is described by the monodromy 
matrix M N (a) — dM N (a) / da . The point a = M(a) 
is elliptic (hyperbolic) if |trMjy| < 2 (|trMjv| > 2). The 
terms elliptic and hyperbolic will be used in the same 
way to describe any trajectory, regardless of whether it 
is closed. 

For chaotic maps, in the large-iV limit the elliptic 
and hyperbolic trajectories are finely intermingled, im- 
plying that P(N) is independent of any smooth scal- 
ing of the phase-space coordinates, and P(N) will de- 
crease rapidly with increasing N. Consider the distribu- 
tion of values of tr Mn(o) for large values of AT. The 
monodromy matrix Mn is a product of elementary mon- 
odromy matrices m(a) = Mi (a) for individual applica- 
tions of the mapping Mat (a) = Y[n=i m ( an ) wnere a n 
is the phase-space point reached from a after n appli- 
cations of the mapping. The typical values of the el- 
ements of the matrix Mn are expected to grow expo- 
nentially, in the sense that the mean of the logarithm 
of some norm of the matrix Mjq should grow linearly 
with time: A = i limjv^oo (log tr [Mjf(ct)Mff (a)]) a 
is termed the Lyapunov exponent. This suggests that 
elliptic trajectories are rare in the long-time limit, and 
that if a universal form exists for the function P(N), it 
might be expected to be exponential, 



P(N) ~ Aexp(-aXN) 



(3) 



for N 3> 1, where A and a are universal constants. 

Dctcrming the fraction of elliptic trajectories is equiva- 
lent to determining the probability that tr Mjv lies in the 
interval [—2, 2]. Because the typical value of trM/v is ex- 
ponentially large, this question relates to the tail of the 
distribution. In the case of a chaotic map the succesive 
positions a n have the characteristics of random numbers, 
and the monodromy matrix Mn may be modelled as a 
product of random matrices. A product of a large num- 
ber of random scalars has a log- normal distribution, and 
it is therefore natural to expect that trMjv will have 
a log-normal distribution. The central limit theorem is 
only applicable sufficiently close to the maximum of the 
probability distribution, and the tails of the distribution 
of a sum of random variables depends upon the distri- 
bution of the variables. Because the distribution of the 
matrices m(a) is non-universal, we therefore expect that 
the tail of the distribution of tr Mn is non-universal. 

Further support for this prediction comes from consid- 
ering moments of Mm {a). Modelling the a n as random 
variables, we have 



(M N (a)) = H(m(a n ))=[(m(a))] 



N 



(4) 



where the averages are taken over the invariant measure 
of the dynamics, which in the Hamiltonian case is just the 
uniform distribution on the phase-space area. Higher mo- 
ments are obtained by averaging outer products of mon- 
odromy matrices: (tr M^(a)Mj\r(a)) for example may 



be expressed in terms of [(m(a) (3 m(a))] (see also M). 
The value of (tr Mjj (a) Mn (&)) is therefore expected to 
grow exponentially with A" 



(tr M%(a)M N (a)) ~ exp(A 2 A0 



(5) 



where A 2 is the largest eigenvalue of (m(a) ® m(a)). A 
similar approach can be used to estimate the growth of 
((tr M^M N ) k ): the result is of the form (5) with A 2 re- 
placed by A 2 fc, the largest eigenvalue of the fc-fold outer 
product (migimig)- • ■0m). The values of these eigenvalues 
depend upon the structure of the elementary monodromy 
matrices m, indicating that the ratios of the eigenvalues 
Afe are non-universal. It is very difficult to reconcile this 
with the hypothesis that P(N) is universal. 

To test these predictions, the fraction of elliptic trajec- 
tories and the Lyapunov exponent were evaluated for the 
mappings defined by Eqs. (1) and (2). In Fig. || we plot 
the fraction P(N) of elliptic trajectories as a function of 
AA^. Fig. clearly shows the non- universality of P{N). 

Next we discuss the nature of the strong fluctuations 
apparent in Fig. [I]. This figure indicates that the den- 
sity of periodic points is highly non-uniform, however 
an exact result due to Hannay and Ozorio de Almeida 
H suggests that the periodic points might have a uni- 
form distribution in phase space. In the following, we 
discuss how this apparent contradiction is resolved, and 
present an argument indicating that the density of peri- 
odic points has a log-normal distribution. 

We base our discussion on the Kac-Ricc jl(|[nj ap- 
proach to calculating the densities of point singularities 
of random functions. The problem is to find the set of 
points where a = M N (a) for large N. The basis of the 
Kac-Rice approach is to estimate the probability 5P of 
finding an example of the singularity in a small ball of 
volume SV centred on an arbitrary test point: this can 
be written SP ~ T> SV, where T> is the density of singu- 
larities. The estimate of SP is usually simplified by the 
fact that the distance to a nearby singularity is given by 
a simple expression. 

Assume that there is a periodic point at a small dis- 
placement 8a from the test point a. Considering the 
condition for periodicity a + 5a = A4(a + 5a), and ex- 
panding the mapping to first order in 5a leads to the 
approximation 



5a = [M -iy\a-M) 



(6) 



(where the index A" designating the AT-th iterate has been 
dropped). For long orbits, it may be assumed that the 
joint distribution function P\M.,M] factorizes, and 



V = P[M] 



M = 

= A' 1 (\det{M - I)\) 



j dM P[M] |det(M — 7)| 

a J 



(7) 



where A is the phase-space area. Similarly, the density of 
periodic points weighted with a smooth function W(a) 
is 
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V w = A' 1 (W(a) \det{M - I)\) 



(8) 



Hannay and Ozorio de Almeida || considered the case 
W(a) = \wj\ 2 , where wj is the weight of the j th peri- 
odic orbit in the Gutzwiller sum. This case is of im- 
portance when estimating quantum two-point correla- 
tion functions semiclassically. The periodic orbit weights 
are of the form Wj = exp(mfij/2) det(M — J) -1 / 2 where 
are integers termed Maslov indices which account for 
the phase changes associated with focusing or reflection. 
This gives W(a) = |det(M — so that in this case 

T>w — ^ L which constitutes the version of the sum rule 
derived in |9fl which is applicable to maps. This simple, 
universal form is a result of the cancellation of the weight 
|det(M — I)\ against its inverse. 

We can now comment on the highly non-uniform dis- 
tribution of periodic points shown in Fig. [I]. It is natural 
to attempt to characterise this by a local density T>(a) 
defined in terms of the number of orbits inside a ball of 
radius, say, e. In order to construct a satisfactory defini- 
tion using this approach, the local density would have to 
converge over a range of values of e, this range becoming 
broader as N increases. The fluctuations in density are 
so wild that it appears to be impossible to define a local 
density in this way. Instead, we define the local density 



V{a) = \6a\ 



(9) 



where 8a is the displacement from the point a to the 
nearest periodic point, \6a\ is the corresponding Eu- 
clidean distance, and d is the phase-space dimension. For 
points sufficiently close to a periodic point, we have al- 
ready seen that \Sa\ ~ |(M - iy 1 (a - M(a))\. The 
fluctuations of this quantity are dominated by those of 
M which are log-normal. We plot the distribution of the 
local density evaluated according to (^|) in Fig. |[ It is 
well described by a log-normal distribution. 

We now discuss the implications of our findings for a 
minor variation upon a classical correlation function of 
periodic-orbit actions first introduced in 



C(T,AS) = ^w J w}6 fl (T-T s ) 
x 5 e (AS - (Sj - Sj,)) . 



(10) 



Here, Sj and Tj are the action and period of the j pe- 
riodic orbit. In systems where Gutzwiller's trace formula 
is exact, and which exhibit universal spectral two-point 
correlations, C(T, AS) must be a function of the action 
difference which is universal within each of Dyson's sym- 
metry classes [||J^] . In general, however, the trace formula 
cannot be exact. It is thus desirable to explore classical 
models for periodic-orbit correlations contained in ([To|). 

A mechanism for such correlations is revealed by 
embedding the Hamiltonian (or mapping) into a one- 
parameter family. We remark that systems for which 
exact trace formulae exist (Riemann's zeroes, motion of a 
surface with constant negative curvature H and quantum 



graphs p2]) do not form one-parameter families exhibit- 
ing bifurcations. The parameter X could be a coupling 
constant such as K in Eq. ([!]) . Varying the parameter X 
produces bifurcations of orbits: immediately after a bi- 
furcation the two orbits have the same actions Sj = Sji , 
and their weights Wj and Wy may be related in a sim- 
ple way (typically they have opposite signs). We will 
estimate the contribution from bifurcations to the corre- 
lation function ([n]) at small values of AS, as this can be 
directly related to the forms of the bifurcations. Assume 
that a bifurcation occurs at Xi as a parameter X is var- 
ied. The action difference in the neighbourhood of the 
bifurcation is of the form 



IAS, 



\AX, 



(11) 



where AXi = X — Xi. The exponent (3 will be deter- 
mined from the type of bifurcation, and the constant Vi 
depends on the particular bifurcation. The weights Wj 
in (10) may be either singular or regular in the neigh- 
bourhood of the bifurcation: we allow for them having 
an algebraic singularity with exponents 7^ , which may be 
different for the two orbits in the vicinity of the bifurca- 
tion, 



Wj (x \AXi\ lj 



(12) 



In the following we discuss a model for evaluating the 
contributions to ( |l0| ) for long periodic orbits, treating 
the Wj in ( |Io| ) as random variables, which are indepen- 
dent except for those pairs of orbits which are related by 
a bifurcation: 



(wj)=0, (wjWj>) 



w 2 5jj> 



(13) 



where in the latter case it is assumed that the orbits j 
and j' are not related by a bifurcation. At a bifurca- 
tion, one or both of the weights may be singular. In this 
case, it is reasonable to model the behaviour of {wjWj>) 
as follows 



J; 2 A -( 7j + 7j ,)| AX |7 3 +7 J 



(14) 



where A characterizes the frequency with which bifurca- 
tions occur. 

Let SP e (AS) be the probability that a periodic orbit 
is connected to another orbit for which the action differ- 
ence is in the interval [AS, AS + e]. According to the 
model above, the only contributions to the correlation 
function come from pairs of periodic orbits related by a 
bifurcation. The contribution from these orbits may be 
estimated as follows 



lfdN\ 2 {AX\ 



(73+7,-') 



6P £ {AS) (15) 



where N(T) is the number of periodic orbits with period 
less than T, and AX is the typical distance to a bifurca- 
tion. From ([ll]) , the bifurcation occurs at a displacement 
in parameter space of the form AX ~ (AS/v) 1 ^ . 
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We model the probability that a periodic orbit does 
not undergo a bifurcation in a distance AX from an ar- 
bitrarily chosen test point as a Poisson distribution: this 
probability is P±(AX) = exp[— A±|AX|] for displace- 
ments to either side of the test point. For a sufficiently 
small separation AX, the probability of finding a bifur- 
cation in a small interval of size SAX on either side of 
X is 



SP ~ A SAX 



(16) 



with A = A + + A_. Combining ( |i~l| ) and ([l6|), one ob- 
tains 6P e (AS) ~ Av^^AS^/^^SAS. Identifying SAS 
with e, and substituting into (^|) produces the following 
result, valid for small |A5|: 



C(T,AS) oc |A5| 



(l+7 3 +7,'-/3)//3 



(17) 



Evaluation of the exponent in ( |17| ) requires information 
about the nature of the bifurcations. As an example, con- 
sider a hyperbolic billiard system with no corners, such 
as the Lorentz gas. In this case inverse bifurcations occur 
when pairs of similar orbits, only one of which bounces 
off a surface, become tangential to that scattering sur- 
face. Geometrical considerations imply that (3 = 2, and 
-fj = for the orbit which does not bounce tangentially 
but 7j-/ = 1 for that which does. In this case the expo- 
nent in (^) is zero, but in other cases the exponent may 
be non-zero. 

We summarise the conclusions on correlations of peri- 
odic orbits as follows. In systems with a smooth Hamil- 
tonian, the statistic P(N) is very closely related to the 
condition for bifurcations (namely that a periodic orbit 
intersects the manifold where trill = ±2). The statistic 
P(N) was shown to be non-universal, and it is implau- 
sible that bifurcations will show universal behaviour in 
such systems. Totally chaotic systems have no elliptic 
trajectories, and bifurcations are associated with singu- 
larities of the Hamiltonian: in this case the arguments 
against universality are no less compelling. We there- 
fore conclude that there is a non-universal component to 
the periodic orbit correlation function. Our result fll7j ) 
also indicates that the contributions to the correlation 
function coming from bifurcations has a singularity if the 
exponent l+7j +7j' — P is non-zero. It is difficult to con- 
ceive of a non-universal and non-analtyic contribution to 
the correlation function which could combine with this 
one to give a universal result. 
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FIG. 1. Shows periodic points of the standard map with 
K = 6 and N = 7 (left) and for the modified cat map K = 2 
and N = 11 (right). In both cases, < x,p < 2n. 
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FIG. 2. Shows P(N) as a function of XN for the stan- 
dard map (K = 15, A ~ log_ftT/2 = 2.015, circles) and for the 
modified cat map (K = 8, A = 1.404, squares). 
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FIG. 3. Shows the distribution of the local density of peri- 
odic orbits X>(«) for the standard map with K = 6 and N = 7 

(a) and for the modified cat map with K = 2 and TV = 11 

(b) . In both cases the distribution is log-normal ( ). 
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